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Abstract: Uniaxially (along Z-axis) anisotropic Heisenberg ferromagnet, in presence of 
time dependent (but uniform over space) magnetic field, is studied by Monte Carlo sim- 
ulation. The time dependent magnetic field was taken as elliptically polarised where the 
resultant field vector rotates in X-Z plane. The system is cooled (in presence of the ellip- 
tically polarised magnetic field) from high temperature. As the temperature decreases, it 
was found that in the low anisotropy limit the system undergoes three successive dynam- 
ical phase transitions. These three dynamic transitions were confirmed by studying the 
temperature variation of dynamic 'specific heat'. The temperature variation of dynamic 
'specific heat' shows three peaks indicating three dynamic transition points. 

PACS numbers: 64.60.-i,05.45.-a,75.60.-d,75.70.-i,75.40.Mg 



Introduction: 

The dynamical behaviours of magnetic model systems, in the presence of time de- 
pendent magnetic field, show interesting physical phenomena [1]. The nonequilibrium 
dynamical phase transition [1], particularly in the kinetic Ising model has drawn much 
interest of researchers in the field of nonequilibrium statistical physics. The dynamic 
transition in the kinetic Ising model in the presence of magnetic field (sinusoidally vary- 
ing in time) was first noticed [2] in the mean field solution of the dynamical equation for 
the average magnetisation. The time averaged magnetisation over a full cycle (of exter- 
nal magnetic field) becomes nonzero at finite values of temperature and field amplitude. 
These values of temperature and field amplitudes depend on the frequency of oscillating 
field. However, the transition there [2] is not perfectly dynamic in nature since it can 
exist for such equation of motion even in static (zero frequency) limit ! This reveals that 
the transition in the zero frequency limit is an artefact of the meanfield approximation 
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which does not consider the nontrivial fluctuations. The occurence of the true dynamic 
transitions for models, incorporating the thermodynamical fluctuations, was later shown 
in several Monte Carlo studies [1]. 

After observing the true dynamic transition in the kinetic Ising model in presence of 
oscillating magnetic field and knowing that it is a nonequilibrium transition, a considerable 
amount of studies were performed [1] to establish that this transition is thermodynamic 
phase transition. The divergences of 'time scale' [3] and the 'dynamic specific heat' [3] 
and the divergence of length scale [4] are two important observations to establish that the 
dynamic transition is a thermodynamic phase transition. 

Although the dynamic transition in kinetic Ising model is an interesting phenomenon 
and a simple example to grasp the various features of nonequilibrium phase transitions, 
it has several limitations. In the Ising model, since the spins can have only two orienta- 
tions (up/down), some interesting features of dynamic transitions (related to the dynamic 
transverse ordering) are missing in this model. The classical vector spin model[5] would 
be the better choice to see such interesting phenomena which are missing in Ising model. 
One of such examples is the 'off-axial' dynamic transitions [6] recently observed in the 
anisotropic Heisenberg ferromagnet. In the 'off-axial' dynamic transition, the dynami- 
cal symmetry along the axis of anisotropy can be broken by applying an oscillating field 
along any perpendicular direction (X-direction say). The dynamic phase transition in 
anisotropic XY spin system in an oscillating magnetic field is recently studied [7] by solv- 
ing Ginzburg-Landau equation. The dynamic phase transition and the dependence of 
its behaviour on the bilinear exchange anisotropy of a classical Heisenberg spin system 
(planar thin ferromagnetic film), is recently studied [9] by Monte Carlo simulation. 

All these studies on the dynamic phase transition, made so far, are related to single 
transition. The dynamic transition occurs at a single value of temperature (for fixed 
values of field amplitude and frequency). No evidence of multiple dynamic transitions 
(for bulk only) is reported so far in the literature in anisotropic Heisenberg ferromagnet 
driven by polarised magnetic field. However, it should be mentioned here, that a very 
recent study [8] of dynamical phase transitions in thin Heisenberg ferromagnetic films 
with bilinear exchange anisotropy, has shown multiple phase transitions for the surface 
and bulk layers of the film at different temperatures. Here, in this paper, the observations 
of multiple (triple) dynamic transitions, in anisotropic Heisenberg ferromagnet driven by 
elliptically polarised magnetic field, are briefly reported which is observed in the uniaxially 
anisotropic Heisenberg ferromagnet (three dimensional) in the presence of a elliptically 
polarised magnetic field studied by Monte Carlo simulations. 

The paper is organised as follows: the model is introduced and Monte Carlo simulation 
technique are described in the next section, third section contains numerical results with 
figures and the paper ends with a summary and few concluding remarks given in the 
fourth section. 

The model and Monte Carlo simulation technique: 

The Hamiltonian of a classical anisotropic (uniaxial and single-site) Heisenberg model 
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[5], with nearest neighbour ferromagnetic interaction in the presence of a magnetic field, 
can be represented as 

H = -jY t §i-§j-DY,{S!Y-h-Y,Su (i) 

<ij> i i 

where Si[Si X , Si y , Si Z ] represents a classical spin vector of magnitude unity (S 2 X + S 2 y + 

Sf z = 1) situated at the i-th lattice site. The classical spin vector Si can be oriented 
in any (unrestricted) direction in the vector spin space. In the above expression of the 
Hamiltonian, the first term represents the nearest neighbour (< ij >) ferromagnetic 
(J > 0) interaction. The factor D in the second term represents the strength of uniaxial 
(z-axis here) anisotropy which is favouring the spin to be aligned along the z-axis. Here, 
it may be noted that for D — 0, the system is in the isotropic Heisenberg limit and 
for D — > oo the system goes to Ising limit. The last term stands for the interaction 
with the externally applied time dependent magnetic field (h[h x ,h y ,h z ]). The magnetic 
field components are sinusoidally oscillating in time i.e., h a = h 0a cosut, where h Qa is 
the amplitude and u> is the angular frequency (cu = 2nf; f is the frequency) of the a- 
th component of the magnetic field. In this present case, the field is taken elliptically 
polarised. The polarised field can be represented as 

h = xh x + yh y + zh z = xh 0x cos(Lut) + zh 0z sm(ujt). (2) 

One can readily check that h x = h 0x cos(ut) and h z = h 0z sm(uit) yield, after the elimina- 
tion of time, 

hi ti 



T~ + JJ~ = 1> (3) 



h 2 h 

"■Ox a 0z 

which shows that the magnetic field lies on X-Z plane and is elliptically polarised (in 
general ho x and ho z are not equal). If h,Q X = ho z = /lo(say), the above equation will take 
the form h 2 x + h 2 z = h\ and the field will be called circularly polarized. The magnetic fields 
and the strength of anisotropy (D) are measured in the unit of J. The model is defined 
on a simple cubic lattice of linear size L with periodic boundary conditions applied in all 
three directions. 

The model, described above, has been studied by Monte Carlo simulation using the 
following algorithm [10]. To obtain the equilibrium spin configuration at a particular 
temperature T, the system is slowly cooled down from a random initial spin configu- 
ration [11]. At any fixed temperature T (measured in the unit of J/K B , where K B is 
the Boltzmann constant), and for the fixed values of ho x , ho z uo and D, a lattice site % 
has been chosen randomly (random updating scheme). Monte Carlo simulations were 
performed using Metropolis algorithm [10] with random updating scheme. The spin-tilt 
trial configuration is generated as follows [11, 6]: two different random numbers r\ and T2 
(uniformly distributed between -1 and +1), are chosen in such a way that R 2 = r\ + r\ 
becomes less that or equal to unity. The set of values of r\ and r2, for which R 2 > 1, 
are rejected. Now, u = y/l — R 2 Si X = 2-ur l5 Si y = 2ur 2 and Si Z = 1 — 2R 2 . In this way, 
the distribution of points of tips of spin vectors on the surface of a unit sphere will be 
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uniform. The acceptance of a trial configuration is determined by Metropolis rate [10]. L 3 
numbers of such updates (at random positions) of spin vectors, defines one Monte Carlo 
step per site (MCSS) and this may be considered as the unit of time in this simulation. 
The linear frequency (/ = uo/2-k) of the time varying magnetic field is taken 0.02 and 
kept constant throughout this simulational study. Thus 50 MCSS are required to obtain 
one complete cycle of the oscillating field. Consequently, 50 MCSS is the time period (r) 
of the oscillating magnetic field. Any macroscopic quantity, such as, any component of 
magnetisation at any instant, is calculated as follows: Starting with an initial random 
spin configuration (high temperature phase), the system is allowed to become stabilised 
(dynamically) upto 4xl0 4 MCSS (i.e. 800 complete cycles of the oscillating field). The 
average value of various physical quantities are calculated from further 4 x 10 4 MCSS (i.e., 
averaged over another 800 cycles). This is important to achieve stable value and it was 
checked carefully that the number of MCSS mentioned above is sufficient to obtain stable 
value of the measurable quantities etc. which can clearly show the dynamic transitions 
points within limited accuracy. But to describe the critical behaviours very precisely (e.g., 
to estimate critical exponent etc.) much longer run is necessary. Here, the total length of 
simulation for one fixed temperature is 8 x 10 4 . The system is slowly cooled down (T has 
been reduced by small interval) to obtain the values of the statistical quantities in the low 
temperature ordered phase. The last spin configuration obtained at previous temperature 
is used as the initial configuration for the new temperature. The CPU time required for 
8 x 10 4 MCSS is approximately 25 Minutes on an Intel Pentium-Ill processor. 



Numerical results: 

The simulational study is done for a simple cubic lattice of linear size L = 20. The 
instantaneous magnetisation components (per lattice site) m x (t) = J2ij%, m y(t) = J2i ~jj 
and m z (t) = are calculated at each time in the presence of magnetic field. The time 

averaged (over a full cycle of the oscillating magnetic field) magnetisation components 
(the dynamic order parameter components) Q x = ^ § m x (t)dt, Q y = \ § m y (t)dt and 
Qz = \ § m z (t)dt are calculated by integrating (over the complete cycle of the oscillating 
field) the instantaneous magnetisation components. The total (vector) dynamic order 
parameter can be expressed as Q = iQ x + jQ y + kQ z . The instantaneous energy e(t) = 
—JJ2<ij> Si- Sj — D Y^i{Si) 2 — h-J2i Si is also calculated. The time averaged instantaneous 
energy is E — ^ § e(t)dt. The rate of change of E with respect to the temperature T is 
defined as dynamic specific heat C{— ^p) [3]. The dynamic specific heat C is calculated 
from energy E, just by calculating the derivative using the three-point central difference 
formula, given below. 

_ dE _ E(T + ST) - E(T - ST) 
C -dT~ 2ST (4) 

For the elliptically polarised (equation 3) magnetic field, where the resultant field lies in 
X-Z plane, the amplitudes of fields are taken as ho x = 0.3 and Hq z = 1.0. The strength of 
uniaxial anisotropy is taken D = 0.2. This value of D is obtained by rigorous searching 
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to have these interesting results and kept constant throughout the study. However, there 
must be variations in transition points depending on the values of D. It is observed that 
for Higher values of D the multiple transition phenomenon disappears. The values of field 
amplitudes and frequency are also obtained by searching. 

The temperature variations of the dynamic order parameter components (Q x , Q y , Q z ) 
are studied and the results are depicted in Fig. 1(a). As the system is cooled down, from a 
high temperature disordered (Q = 0) phase, it was observed that, first the system under- 
goes a transition from dynamically disordered (Q = 0) to a dynamically Y-ordered (only 
Q y 7^ 0) phase. This may be called as the first phase (Pi) and the transition temperature 
is T cl . This phase can be characterised as P^ {Q x = 0, Q y 7^ 0, Q z = 0). Here, the 
resultant vector of elliptically polarised magnetic field lies in x-z plane and the dynamic 
ordering occurs along y-direction. So, this is clearly an off-axial transition [6]. In the case 
of this type of off-axial transition the dynamical symmetry (in any direction; y-direction 
here) is broken by the application of the magnetic field in the perpendicular direction 
(lies in the x-z plane here). As the system cools down, it retains this particular dynami- 
cally ordered phase (Pi) over a considerable range of temperatures. As the temperature 
decreases further, a second transion was observed. Here, the system becomes dynami- 
cally ordered both in X- and Z-directions at the cost of Y-ordering. In this new dynamic 
phase, P 2 : (Q x 7^ 0, Q y = 0, Q z ^ 0). In this phase the dynamical ordering is planar (lies 
on x-z plane). The ordering occurs in the same plane on which the field vector lies. This 
transition is axial [6]. This phase may be called the second phase (P2) and the transition 
(from first phase to the second phase) temperature is T c2 . As the temperature decreases 
further, the X- and Z-ordering increases. At some lower temperature, a third transition 
was observed, from where the X-ordering starts to decrease and only Z-ordering starts 
to increase quite rapidly. This third phase can be designated as P3: (Q x 7^ 0, Q y = 0, 
Q z 7^ 0). Although the characterisation of P2 and P3, in terms of the values of dynamic 
order parameter components, looks similar there exists an important difference between 
these two phases. In the phase P2, both Q x and Q z increases as the tempereture decreases 
but in the phase P 3 , Q x decreases as the temperature decreases (see Fig. 1(a)). So these 
two phases P2 and P3 distinctly differs from each other. In this phase the dynamical 
ordering is also axial (along Z-axis or anisotropy axis). The system continues to increase 
the dynamical Z-ordering as the temperature decreases further. The low temperature 
phase is only dynamically Z-ordered. That means the systems orders dynamically (only 
Q z 7^ 0) along the Z-direction (direction of anisotropy) only at very low temperatures. 
Zero temperature dynamic phase (for such polarised field) can be characterised as Q x = 0, 
Q y — and Q z = 1.0. 

To detect the dynamic transitions and to find the transition temperatures the tem- 
perature variation of the energy E is plotted in Fig. 1(b). From this figure it is clear 
that there are three dynamic transitions occur in this case. The transition points are 
the inflection points in E — T curve. The temperature derivative of the energy E is the 
dynamic specific heat C. The temperature variation of C is shown in Fig. 1(c). The three 
dynamic transitions are very clearly shown by three peaks of the specific heat plotted 
against the temperature T. From this figure the transition temperatures are calculated 
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(from the peak positions of C — T curve). First transition (right peak) occurs around 
T c i = 1.22, the second transition (middle peak) occurs at T c2 = 0.94 and the third (left 
peak) transition occurs around T c3 = 0.86. 

This study was further extended for other values of h Qx keeping other parameters fixed. 
It was found that this three transitions senario disappears for higher values of Hq x . For 
example, for h$ x = 0.9, the second phase P<i disappears. In this case, the C — T curve 
shows two peaks. It was also observed that for h 0x = 0.2, h 0z = 0.2 (keeping all other 
parameter fixed) the system shows single transition and only dynamically orders along 
Z-direction. 

To detect the dynamic transition points an alternative method may be to study the 
temperature variation of fluctuation of dynamic order parameter x{Q)[ = L 3 (< Q 2 > — < 
Q > 2 )]. However, we do not have sufficient amount of precise data to study this. 

Summary: 

The uniaxially (Z-direction) anisotropic Heisenberg ferromagnet in presence of time 
dependent (but uniform over space) magnetic field is studied by Monte Carlo simulation 
using Metropolis dynamics. The dynamic transition in uniaxially anisotropic Heisenberg 
ferromagnet is already studied by MC simulation. In that case the time dependent mag- 
netic field was sinusoidal and the axial and off-axial dynamic transition was reported [6] 
earlier. 

In the present study, the external time dependent magnetic field was taken elliptically 
polarised where the resultant field vector rotates on X-Z plane. For the lower values of 
anisotropy and a specific range of the values of field amplitudes the system undergoes 
multiple dynamic phase transitions. Here, three distinct phases are identified. In this 
paper, this observation is just briefly reported. This multiple dynamic phase transition in 
anisotropic Heisenberg ferromagnet in presence of elliptically polarised field, is observed 
here by Monte Carlo simulation. An alternative method, to check this phenomenon, 
may be to use Landau-Lifshitz-Gilbert equation of motion [12] with Langevin dynamics. 
Another important thing should be mentioned here regarding the possible explanation 
of multiple dynamic phase transitions (axial and off-axial transitions) observed in the 
anisotropic Heisenberg model. One possible reason may be the coherrent rotation of 
spins. Where the dynamic phase transition in the Ising model can be explained simply 
by spin reversal and nucleation [13]. But to establish the responsible mechanism behind 
the multiple dynamic phase transitions, details investigations are required. 

The variations of the dynamic phase boundaries with frequency and the strength of 
anisotropy is quite interesting to be studied. This study also indicates that the system will 
show a very rich phase diagram with multicritical behaviour. The finite size analysis is 
also necessary in order to distinguish the crossover effects from the true phase transitions. 
This requires huge computational task which will take much time. This work is in progress 
and the details will be reported later. 
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Fig. 1(a) 

Fig. 1(a): The temperature variations of the components of dynamic order parameters. 
Different symbols represent different components. Q x (O), Q y ( Q ) and Q z (»)• This 
diagram is for D = 0.2 and for elliptically polarised field where ho x = 0.3 and Hq z = 1.0. 
The size of the errorbars of Q x , Q y and Q z close to the transition points is of the order 
of 0.02 and that at low temperature (e.g., below T = 0.5) is around 0.003. 
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Fig.l(b) 

Fig. 1(b): The temperature variation of the dynamic energy (E) for D = 0.2, h 0x = 0.3 
and h 0z = 1.0. The vertical arrows represent the transition points. 
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